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HYPERSONIC LAMINAR VISCOUS INTERACTION EFFECTS 
ON THE AERODYNAMICS OF TWO-DIMENSIONAL WEDGE 
AND TRIANGULAR PLANFORM WINGS 

By Mitchel H. Bertram 
Langley Research Center 

SUMMARY 

Solutions have been obtained of the effect of viscous interaction not only on the pres- 
sure and skin friction but also on the aerodynamic forces of a sharp-leading-edge flat sur- 
face. The analysis assumes the applicability of laminar-hypersonic-local-similarity- 
boundary-layer theory in a perfect gas with Prandtl number unity and a constant ratio of 
specific heats. The results are presented in the form of correlation graphs which cover 
a large range of plate incidence and viscous interaction parameter and allow the effects 
of viscous interaction to be readily determined. Computations are shown for ratios of 
specific heats of 7/5 and 5/3 and, for the aerodynamic forces, both the two-dimensional 
wing and the triangular planform wing are considered. The results indicate that the effect 
of viscous interaction can reduce the lift -drag ratio significantly. 

INTRODUCTION 

Utilizing laminar-hypersonic-local-similarity-boundary-layer theory yields rela- 
tively simple solutions in correlation form and allows the pressure, skin friction, and 
aerodynamic coefficients for sharp plates or wedges to be obtained. The analysis which 
follows is based on the work presented in references 1 and 2 and applies for Prandtl num- 
ber unity in a perfect gas with a constant ratio of specific heats. In earlier papers by the 
present author and his coworkers (refs. 2 and 3), the effect of angle of attack on the 
boundary-layer-displacement-induced pressures was taken into account by assuming the 
boundary layer grows in the inviscid flow behind the plane shock produced by the plate. 
With this approach it was found that the viscous effects on lift and drag essentially can- 
celed and the lift-drag ratio was relatively unaffected by displacement effects. White, 
however, has shown the more nearly correct approach is to consider the flow to be dis- 
placed by both the plate inclination and the boundary layer. (See ref. 4.) In this case the 
foregoing effect on lift-drag ratio cannot be assumed to apply and the calculation of the 
effect on surface pressure, skin friction, and aerodynamic forces on a flat plate or wedge 
is the subject of this paper. 



SYMBOLS 


a constant in expression for (see eq. (4)) 

a,b,c coefficients in approximation for fw (see eqs. (33)) 

c local streamwise chord of triangular planform wing 

c r root chord of plate with triangular planform 

C coefficient in linear formula for viscosity (see eqs. (A4) to (A6)) 

C ^ chord force coefficient due to pressure 

Cjj drag coefficient 

Cf local skin-friction coefficient 

C F average skin-friction coefficient 

C L lift coefficient 

C m moment coefficient about leading edge of flat plate or apex of delta planform 

wing (M/qSL), C m>A = M/qSc r 

C N normal-force coefficient 

C N a H dC N/ dQ! 

f^' shear stress function 

G function of wall temperature and specific heat ratio in laminar-boundary- 

layer growth equation (see appendix A) 

K total effective flow deflection angle in hypersonic similarity form, K Q + 

K q angle of plate surface relative to free stream in hypersonic similarity form 

positive for windward facing surface, negative for leeward facing 
surface, 
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coefficient in laminar-hypersonic-local-similarity-boundary-layer theory 
local skin-friction equation (see eq. (31)) 

coefficient in hypersonic-similarity-theory-laminar-boundary-layer growth 
equation (see eq. (4)) 

angle of attack of wing in hypersonic similarity form, 

local inclination of boundary layer with respect to plate surface in hypersonic 
similarity form, s M^dfi/dx 

length of plate 

lift-drag ratio 

moment 

free-stream Mach number 

short notation for value of (P - Po)A as X — ®° 

Prandtl number 

exponent in power-law variation of surface pressure with surface distance 

ratio of local surface pressure to free-stream static pressure 

ratio of inviscid surface pressure to free-stream static pressure 

ratio of average pressure on one side of plate of length L to free-stream 
static pressure 

ratio of average pressure on one side of plate of triangular planform with 
root chord c r to free-stream static pressure 

dynamic pressure 


Reynolds, number 



Xtl 


s 


planform area of wing 


T temperature (absolute) 

x streamwise distance on plate measured from leading edge 

x cp x-location of center of pressure 

a angle of attack of wing 

j3 pressure gradient parameter in transformed plane of laminar-boundary- 

layer similarity theory 

y ratio of specific heats 

6 boundary-layer thickness 

6 angle of inclination of surface relative to free-stream flow 

X viscous interaction parameter, Gyj 2 

/i dynamic viscosity 

p fluid density 

r half-angle of wedge section airfoil 

viscous interaction parameter, \[c/fR ~ L 

Subscripts: 

c based on root chord 

r 

l lower surface of wing 

L based on length L 

o two-dimensional case without viscous interaction effects 
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r adiabatic wall 

t stagnation 

T total (including both surfaces of wing) 

u upper surface of wing 

w wall 

x based on distance x 

°° based on undisturbed free-stream conditions 

A triangular planform 

lim limiting value 

max maximum 


ANALYSIS 


Pressure Distribution 


Basic to this analysis is a solution for the pressure distribution on a plate at arbi- 
trary angle of attack. The inclination of the plate surface in hypersonic similarity form 
is M^e^Ko) and the local inclination of the boundary layer with respect to the plate sur- 
face is M oo 6^(=K 6 ). Thus, since the effective turning angle of the fluid is assumed to be 
the sum of these two angles, in hypersonic similarity form K = Kq + Kg. The relation 
Kg is given by equation (7) of reference 2 as 


K 


5 y'p 2 


+ ±&) 

2P dx) 


This relation can be more generally written (so that calculations are independent of 
wall temperature ratio) as 


K* = K, -= 
6 4 


-4(i + 


_X_ dP 
2P d\ 


( 1 ) 


where X is the notation of reference 4 for Gy/2 and G is a simple function of wall 
temperature and specific heat ratio to be given later. The solution of equation (1) requires 
a relationship between P and K and this relationship is assumed to be given by the 
hypersonic similarity shock and expansion equations (tangent wedge theory): 
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(K = 0; P = 1) (2) 



o 

and for K = — — it is assumed that P = 0. Based on the values for the boundary- 

y - 1 

layer-growth coefficient given in figure 1 of reference 2 (modified from those given 

in ref. 1), is approximated by 


= 1 + an (4) 

where n is the local exponent in a power law fit to the surface pressure variation given 
by n = -(XdP/dX)/2P and a is assumed to be independent of wall temperature ratio 
with values as follows: for y = 7/5, a = -0.19; for y = 5/3, a = 0. 

To obtain the pressure distribution the system of equations (1), (2) or (3), and (4) 
has to be solved. Such solutions were previously obtained by iteration in reference 2 for 
certain simple cases but in the present case were made on a high-speed digital computer 
and some details of this solution are given in appendix A. 

An exact closed form solution is obtainable, however, in the limit as X — 0 
where dP/dX = (P - P Q )/X = (dP/dK) dK § y/dX. Then from equation (1), dK^/dX - l//P^ 
For a windward-facing plane surface, equation (2) is utilized and the result is 



P - Pq\ y 

x /X=0 





+ 



(Kq^O) 


( 5 ) 


where P 0 is obtained from equation (2) with K = Kq. For Kq = 0, equa tion (5) reduces 
to (P - P 0 )/X = y and for Kp = °°, the result is (P - P Q )/X = \fi y(y + 1). 

In the case of a plane leeward surface, with equation (3) the solution is 



where Kp is taken as positive for windward surfaces and negative for leeward surfaces. 
Accurate values from equations (5) and (6) and the inviscid pressure ratios from 
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equations (2) and (3) are given in table I. A number of approximate solutions to the pre- 
ceding system of equations are possible and some of these are discussed in appendix B. 

The results from the high-speed digital computer solution of the system of equa- 
tions (1) to (4) including the limiting case obtained from equations (5) and (6) are given in 
figure 1 for a large range of plate incidence and values of y of 7/5 and 5/3. The form 
of pressure parameters is that suggested by the asymptotic solutions and allows accurate 
values of pressure ratio to be obtained from a plot which covers an extremely large range 
of X. 

Values of the exponent n obtained in conjunction with the solution for pressure are 
shown in figure 2. Knowledge of this exponent is necessary for proper evaluation of skin 
friction and heat transfer. (See ref. 1.) Clearly, the values of n can differ markedly 
from the value of -1/2 given by strong interaction theory even for relatively large values 
of X. For large angles of incidence of the leeward surface (negative Kq), n tends to 
have large negative values which in some cases exceed the limits of applicability of the 
hyper sonic -boundary -layer similarity theory. This problem is discussed in a later sec- 
tion where the skin friction is evaluated. 

Pressure Forces 

If it is assumed that the induced pressure distribution on the plate is known from the 
analysis in the preceding section, the effect of the boundary layer on the aerodynamic 
pressure forces and skin friction may be determined. In this section only the aerody- 
namic pressure forces are treated. 

Average pres sure ratio on a flat plate .- If the average pressure over the plate is 
known, the normal and chordwise pressure forces are readily determined. The average 
pressure ratio on a two-dimensional plate is defined as 


P = 



P dx 


(7) 


where x is streamwise distance on the plate measured from the leading edge and L is 
the value of x at the trailing edge of the plate. In terms of the dimensionless boundary- 
layer interaction parameters, equation (7) becomes 


P » 2X L 2 



( 8 ) 
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This statement of P has drawbacks because P (and P) — °° as X — 0. There is, 
however, a convenient way to avoid this problem. Equation (8) can be written in the fol- 
lowing form: 


P-Pq 

X L 




(9) 


Values of (P - P G )/X vary over a relatively small range and never exceed 0(1) as shown 
in figure 1. For very small values of X, the values of (P - P 0 )/X are close to those 
given for the X = 0 solution in the previous section as equations (5) and (6). For X — °°, 
(P - P 0 )/X approaches a constant which is actually the strong interaction solution; thus, 
some further simplification is allowed for machine solution. Equation (9) is now written 
as 



as X — 0, (P - P 0 )/x L — 2(P - (See eqs. (5) and (6).) The general strong inter- 

action solution given in appendix B gives for the first term in equation (10) 



where (K4) n __i/2 * s obtained f rom equation (4). The Xj values are arbitrarily large 
values of X for which P » P Q . The numerical machine solution of equation (10) (with 
eq. (11)) utilized Xj = 10^. This solution is shown in figure 3. Here the average pres- 
sure parameter is shown as a function of X for various surface inclinations. 

The triangular planform is also of interest as a lifting device. In this case the 
assumption is made (as in ref. 5) that each filament in the stream and boundary layer 
over the wing remains unaffected by adjacent filaments and thus that any section of this 
wing may be treated as a section of an unswept plate with the same local chord. The 
average pressure ratio on the triangular plate is thus given by a spanwise integration of 
the average pressures from the previous two-dimensional solution and is expressed in 
terms of the root chord c r of the triangular planform 

= o r c r — 

P =~ \ Px dx (12) 

c r " J 0 
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where x is measured along the root chord from the apex of the wing. By following the 
same type of development as for the two-dimensional case, equation (12) may be written as 



Po 



(13) 


where X r is the value of X based on the root chord of the triangular planform plate 
and j(P - P 0 )A) x ^ o is given, as before, by equation (11) and equation (13) is solved 
numerically. As X — 0, ^P - Po)/^c r 8/3 (P - Po)/^j/ (See eqs. (5) and (6).) This 
solution giving the average pressure parameter on the triangular planform plate as a 
function of the viscous interaction parameter X for a large range of similarity incidence 
angle is presented in figure 4. 

Normal force .- These average pressures allow the determination of normal- and 
chordwise -force coefficients as follows: 

For normal force by definition; 

two-dimensional planform with viscous interaction: 

C N = 1 ~ Pu) 2 / yM °° 2 

triangular planform with viscous interaction: 

c n A ' ' P u) 2 A m » 2 

inviscid: 

C N,o = ( P o,Z " P o,u) 2 / yM °° 2 

where the subscript l refers to the lower surface of the wing and the subscript u to 
the upper surface of the wing. 

These equations may be rearranged into the following convenient forms: 
for the two-dimensional case, 


(14a) 


(14b) 


(14c) 


c c (Izlo) Jlzh o) 

C N ~ _ \ A. ii \ X /u 

x L C N,o P o,Z “ P o,u 


(15) 
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mill 


and for the triangular planform, 

C XT a - C x . 

N,A N,o 

X c r C Nj0 

The values of the pressure parameter for use in equations (15) and (16) are obtained from 
equations (10) and (13) (inviscid pressure ratios from equations (2) and (3) with K = Kq). 
In addition, for X — 0, 

C N,A ~ C N,o 4 C N " C N,o 

^c r ^N,o ^ ^-L^NjO 

An example of the normal-force-coefficient parameters determined from equations 
(15) and (16) for plates of zero thickness is shown in figures 5 and 6 where the parameter 
is shown as function of similarity angle of attack. To aid in interpolation, the result at 
zero angle of attack is shown in figure 7 although X has been removed from the normal - 
force parameter in order to show more directly the effect of viscous interaction on the 
initial normal -force -curve slope. 

Chord force. - For the chordwise force of a wedge-section airfoil by definition, with 
base pressure zero, 

two-dimensional planform with viscous interaction: 



C A = ( 5 i * Pu^M- 2 
triangular planform with viscous interaction: 

C A,A - (P; * Pu^M- 2 

inviscid: 



(17a) 


(17b) 


(17c) 


where r is the half-angle of the wedge section. 

These equations may also be rearranged into the following convenient form: 
for the two-dimensional case, 
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C A " C A,o \ X )l 


P^Po 


^L^AjO 


and for the triangular planform, 


C A,A ~ C A,o 
x c r c A,o 


p o,Z + Po,u 


^Po + P^Po 


P o,Z + P o,u 


The values of the pressure parameter for use in equations (18) and (19) are obtained from 
equations (10) and (13) (inviscid pressure ratios from equations (2) and (3) with K = Kq). 
In addition for X — 0, Kq — 0, 


C.a-C* C, - C* 

•A, A A,o _ g A A,o ___ ^ 

Xc r C A,o X L C A,o 


Moment coefficient and center of pressure .- In coefficient form the moment about the 
leading edge of a two-dimensional flat plate is 


( p * ■ Pu )* dx 

or in terms of the boundary-layer-interaction parameters 

m y L Js \ 5 


y u J x T x5 


For the inviscid case, 


C m,o M oc 


2 P o.l ~ P o,u 


and combining equations (21) and (22) yields 


AC M 2 = (C - C )M 2 = - X T 4 
m 00 ' m m,cr 00 y L 


r m- 

Jx T x 4 


P - P 0 N 

X , 


For present purposes the most useful solution of equation (23) is obtained by separating the 
contribution of the upper and lower surfaces 


AC m M - 2 = (AC m M TO 2 ) 


lAC m M ro 
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where 



( 25 ) 


As in the solution to obtain average pressure, further simplification for numerical solu- 
tion is possible by utilizing the strong interaction solution to obtain the first part of the 
integral extending from X = °° to X = Xj. (See, for example, eqs. (9) and (10).) In the 
present case, 


! AC M ‘ 

m 00 


X L 


l ,u 



+ 




(26) 


with [(P - P 0 )/x) A _oo given by equation (11). In the limit as X — 0, the solution becomes 


^ C m M ~ 


l ,u 
X=0 


P - P, 


x=o 


(26a) 


where [(P - P o )A]x = 0 i s obtained from equation (5) or (6). 


The individual contributions of the windward and leeward surfaces of a flat plate to 
the moment coefficient are given in figure 8 as function of X for a wide range of plate 
inclination. If the indication of the asymptotic results in equations (26) and (26a) are uti- 
lized, however, one finds that a good approximation to the moment coefficient contribution 
at all X and Kq is 


AC M 
m 


2 

00 


X L 



P - P 
*0 

X 


(26b) 


The center of pressure (measured from the leading edge) is given by 


Eq. (26) or (26b) 


Eq. (22) 


\ 


X C n CM' 
1X1 00 


_( 


2^ . „ I ,, 2 


L c n m J 


AC M 2 V - (AC M 2 ) +C M 
m °q jl ^ m /u m,o «> 


'N,o m « 


C N ~ Cn,q 
x L c N,o 


(27) 


Eq. (14c) 


Eq. (15) 
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The effect of viscous interaction on the center-of-pressure location (independent of y) on 
a two-dimensional flat plate is shown in the following sketch. In general, as viscous inter- 
action effects increase (increasing X), the center of pressure moves toward the leading 
edge from mid chord. 


.5 


cp 

L 


. 4 ' 

10 


-3 



As X ■* 0 °, x C p^L — 1/3 and this result would seem to be at variance with the 
trend of the machine solution at large X shown here. However, this apparent discrep- 
ancy is not believed to be important as the theory is generally not valid at the largest X 
and the machine solution is believed to be accurate in the valid range of X. (See a later 
section on "Limitations of Theory.") 

In coefficient form, the moment about the apex of a delta planform wing is 


C 


m,A M ~ 


JL C Cr 

r Cr J 0 




(28) 


and for the inviscid case 

c m,o,A M oo 2 = ( P o ,1 " P o,u) ( 29 ) 

where c is the local chord of the wing measured streamwise, c r is the root chord, and 
x C p^c is the nondimensional streamwise center of pressure of the local chord element 
based on the two-dimensional calculation (eq. (27)). In terms of the viscous interaction 
parameter and with the use of the same procedure as for the two-dimensional case, com- 
bining equations (28) and (29) yields 
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which is equivalent to equation (23) for the two-dimensional case. A running integration 
of equation (30) could not be performed; therefore, because of machine time and complexity 
considerations, equation (30) was not included in the solutions presented in this paper. 

Skin Friction 

Local skin-friction coefficient from hypersonic laminar -boundary-layer theory as 
given in reference 1 is 


C f = 0.664 K 




PC 


R x ,oo 


(31) 


where C is the coefficient in the linear law for viscosity and Kj is the coefficient 
which accounts for the effect of pressure gradient on the local skin friction. Also, Kj 
is a function of wall temperature ratio and the value of the exponent n in a power-law 
variation of wall pressure with distance from the leading edge. In the present case local 
similarity is assumed; thus local values of n are used. (See statements after eq. (4) 
and fig. 2.) Values for = f(n,T w /T 0 ) are given in reference 1; however, for machine 
computation it was convenient to determine Kj as follows: 


K = w 

1 0.4696 


vT 


+ n 


(32) 


(eq. (7) of ref. 1) where fw was approximated by the following equations: 


— — — = 1 + a/3 + b/3 2 + c0 4 
0.4696 


a = 0.627 + 1.8301-=-— 
A t/ 


-w 


= - 0.312 + 0.662^~j + 0.038^~| 

T w /T t 


c = 0.0600e 

and /? and n are related by 


5/2 


(33) 


0 = 


1 - y n 

y n + 1 


(34) 


(eq. (1) of ref. 1). Equations (33) represent an error of less than 1 percent compared with 
exact values for the range 0 S j3 = 1.2; and were originally used for all the calculations. 
However, it was found, after calculations were done, that for the larger negative values of 
angle of incidence Kq, 0 could exceed 1.2 by large amounts and thus could have 
extremely large errors. It is believed that local similarity is not valid for these cases; 
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however, in these cases the pressure is very low and it should be necessary only to insure 
that Kj is of a reasonable magnitude. Thus equations (33) were used for 0=0 = 1, and 
for 0 > 1, this equation was extrapolated linearly with n utilizing the slope obtained 
from equations (33) with equation (34) at 0 = 1 or 


where 


KW - taW 




ill = i J&Y - 

n/fel vV y - 1 L 


a + 2b + 4c - 


2(2 r - l) 


(1 + a + b + c) 


0 


(35) 


and a, b, and c are as given in equations (33). 

Average skin friction, two dimensional .- For one side of a two-dimensional plate, 
the average skin friction is defined as (and with equation (31)) 


C 


F " 



dx = °- 664 ^~ f R °°> L 
Roo,l 



dR 


x,' 


(36) 


If equation (36) is put into terms of the viscous interaction parameter and normalized with 

respect to the skin friction on a flat plate at zero incidence without viscous interaction 

C P , there is obtained 
r ,o 



(37) 


Now P -► °o as X — °° but this problem may be circumvented for purposes of numerical 
integration in a manner analogous to that utilized in solving for P and P in previous 
sections. In other words, a closed form solution will be obtained for X > Xj and a term 
determined by numerical integration will be added for X < Xp Thus, equation (37) may be 
written as 


Cp 

c F,o 



Kjv/P 

Jx L ^r 


dX 


(38) 


If the flat-plate equation (eq. (11)) is used for the value of P in the left-hand integral, 
the solution for equation (38) is with [(P - P 0 )/x]x— 00 = m 
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'jPp + mXj + \[p~q 
^P Q + mXj - /Pi 



(39a) 


For P Q = 0 or very small, equation (39a) is inadequate, but for P 0 « mX, a 
solution is readily determined: 




\2~ 


dX 


(39b) 


In the machine calculation, equation (39a) was used when P Q ^0.001 and equation (39b), 
when P 0 < 0.001. Actually, there is little difference in the result if equation (39b) is 
used for all the calculations regardless of the value of P Q . For example, if Xj is 100 
times greater than P D , the left-hand term of equation (39b) is only 0.1 percent different 
from the left-hand term of equation (39a); if Xj is only 10 times greater than P Q , the 
difference in the left-hand term of the two equations is still less than 1 percent. 

Results from the high-speed digital computer solution of equations (39) are shown in 
figures 9 and 10. Figure 9 presents results for the plate surface at zero incidence and 
figure 10 shows the results for windward and leeward surfaces taken individually. Values 
of \fP^ useful in determining the asymptotic values of the skin-friction ratio (X — 0) are 
given in table I. 

Average skin friction - t riangular planform plate.- If the development in appendix C 
of reference 5 is followed, the average skin friction for one side of a plate of triangular 
planform may be written, the viscous interaction parameter being used, as 

C F,A = 4x c r 4 j! C F x5 ^ 

Ac r A 

where, as before, X is based on the streamwise local chord and X„ is a value of X 

< -r 

based on the root chord. A general solution is readily obtained in the following form: 


C F,A 

^F,o,c r 


= 4Xc t 


i; 


'F,o 


dX 

X4 


(41) 


where _ _ is the average skin-friction coefficient on a two-dimensional plate at 
r ,o,c r 

zero incidence, whose length is the same as the root chord of the delta wing, without vis- 
cous interaction. Values of Cp/Cp Q are obtained from the previous solution of equa- 
tions (39). For numerical integration in this form, large values of X are a problem; as 
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before, equation (41) is rewritten as follows: 


C F,A = 4X 3 / f °° _^F_ dX + C X i _^F_ dx\ 
c F,o,c r Cr C F,o *- 4 " x c r C F,o 


(42) 


As for the two-dimensional case (eqs. (39)), the first part of the integral can be evaluated 
in closed form. For simplicity sake, only the small P Q solution was made (equivalent 
to the solution given as equation (39b)). As noted in the remarks after equations (39), 
this simplification actually affects the accuracy of the result very little for the range of 
K q considered. Thus from equation (39b), since ^yJ ^¥ ,o as 

X equation (42) may be written as 




c F,o,c r 



4 ( K l)n=-l/2 fm + ^ C F dX 
5 Xi 2 V x i ~X Cr C F>0 X 4 


(43) 


as C FAl C F,o,c r - {i/3) ^' 

Results from equation (43) by the high-speed digital computer are shown in figures 9 
and 11. Figure 9 presents results for the triangular planform plate surface at zero inci- 
dence and figure 11 shows the results for windward and leeward surfaces taken 
individually. 

Total skin friction of wedge-section airfoil wing .- Summing up the contribution from 
each of two surfaces of the wing yields for the two-dimensional case 


C F,T 

C F,o,T 




(44) 


where C F Q T = 2C p Q and C F = 1.328 V^/V r l ,°° and values of C F /c F>Q are 
obtained from equations (39). 

Similarly, for the triangular planform wing, 


c f,a,t 

C F,o,c r ,T 


1 

2 


V F,o,c r / l 


'F,A 


'F,°,Cr. 


u 


(45) 


where C F}0 ,c r ,T = 2C F,o,c r and C F,o,c r = 1-328 
^Fja/^FjOjCp are obtained from equation (43). 


and values of 
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The Lift-Drag Ratio 


The lift and drag coefficients for a two-dimensional wedge-section airfoil of half- 
angle r are as follows: 


C T = cos a- C A + C 


'L " N 


A F,T ,, 2 

r M co > 


4r \ . 
sin a 


C~ = C XT sin a + C A + C 


4t 


'D ~ N 


A F,T 2 


icos a 


(46) 


where C^, the first term in parentheses, is from the previous analysis in which the base 
pressure is taken as zero and the last term in parentheses corrects the base pressure to 
that of the undisturbed free stream. Expressing equations (46) as a ratio yields the lift- 
drag ratio 


2 tan a 


L C N'( C A tC F,T' rM> 

D 0 4t 

Cjj tan a + + C p, T 2 


(47) 


in which the parts are readily separable into terms with and without viscous interaction 
as follows: 


C N " C N,o 


1 + Xt [ C N~ C N>°] 


where equations (14) and (15) are used and 


C A = C A,o 


1 + Xt 


x L C N,o J] 


( C A - ^A.oi 


L \ x L c A,o 

where equations (17) and (18) are used, and C F T is obtained from equation (44) with 

equations (39) and Cp,o,T = 2.656 / U j {R^L- Equations (46) and (47) are adaptable to a 
plate of triangular planform by an appropriate change of subscripts; change C N to C N A , 

C'A ^A A’ T ^F,A,T’ X L ^ x Cj>’ an< ^ ^°°,L 

EFFECT OF VISCOUS INTERACTION ON THE LIFT-DRAG RATIO 
OF WINGS WITH ZERO THICKNESS 


The solutions so far presented can generally be applied to wedge-section airfoil 
wings with arbitrary angle although certain examples have been computed which 
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specifically apply to plates of zero thickness. (See figs. 5 to 7 and sketch on page 13.) 
When the lift-drag ratio is computed, individual cases are readily done with the use of 
these graphs; however, the computation of a large number of cases can increase the prob- 
lem of preparation and presentation to unmanageable proportions. Thus, in the present 
case only the results for the zero thickness wing are considered. 

The calculations were done for values of y of 7/5 and 5/3, the two-dimensional 
and triangular planform, and wall temperature to stagnation temperature ratios in the 
range 0 to 1. Required values of G and C were calculated as shown in appendix A. 

The values of total average skin friction at zero angle of attack without viscous interaction 
were found to be as shown in figure 12. 

An example of the sort of curves from which (L/D) max was obtained is shown in 
figure 13(a) where the lift -drag ratio is shown as a function of lift coefficient for the par- 
ticular case of a two-dimensional wing in helium flow at a Mach number of 20 with adia- 
batic wall conditions. In this case the effect of viscous interaction is marked even at high 
Reynolds numbers. This same case is shown in figure 13(b) but the lift coefficient is 
shown as a function of the hypersonic similarity angle of attack. At very small angles of 
attack where and are essentially linear with a, accurate values of can 

be obtained from the equation 


_3 

' L ,a-0 M 


a 


f dC N C N 


, d “ C N,oJ^o \ C F,o/^ 0 F ’°’ T 


(48) 


where, for a flat plate consistent with the hypersonic approximations so far used 
^dC^da^Q = 4/M^ and values of c n/ c n 0 and C f/^F o at 01 = 0 are obtained 

from figures 7 and 9. The value of C„ T can be obtained from figure 12 or any other 
laminar theory. Note that in the example presented in figure 13 (b), the curve for 
without viscous interaction is linear up to about K a = 0.4 and for the cases with viscous 
interaction the extent of linearity increases as Reynolds number decreases until for X = 5 
the curve is essentially linear up to about K a = 1. 

The overall results for the effect of viscous interaction on (L/D) max are shown in 
figure 14 where the ratio of (L/D) max with viscous interaction to the value of (L/D) max 
without viscous interaction is shown as a function of the viscous interaction parameter X. 
One readily sees that there is a significant effect of wall temperature on the importance of 
viscous interaction and that viscous interaction can reduce (L/D) max significantly. The 
effect of temperature is ascribed to pressure-gradient effect on the skin friction since in 
the present context there is no effect of wall temperature on normal force when X is 
utilized as the variable. At the highest Reynolds numbers (lower values of X), the effect 
of viscous interaction, is more or less independent of Mach number but is Mach number 
dependent at lower Reynolds numbers. On the other hand, the changes in (L/D) max ratio 
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due to the change in planform from two dimensional to triangular in shape or a change in 
y from 7/5 to 5/3 are secondary. A plot such as figure 14 is informative but is restric- 
ted in usefulness. Linear theory suggests another type of presentation since according to 
this theory 

(L/D) max 


A correlation of the present results (C^ = 0) according to the parameters suggested 
by equation (49) is shown in figure 15. The curves representing the (L/D) max param- 
eter were found to be independent of wall temperature; however, they were a function of 
Mach number. For the cases with viscous interaction the results correlate both for the 
two values of y considered and for the two types of planform. For the cases without 
viscous interaction, the results correlate for the two types of planform considered but 
differ somewhat for the two values of y. 

From figure 15, values of (L/D) max for a wide range of conditions can be found. 
To obtain the required values of c Nq,/ C F at 01 = 0 with viscous interaction, figures 7 
and 9 could be used but it has been found expeditious to combine these two and the result 
is presented in figure 16. Thus with figure 16 only a knowledge of Cjj > 0 (in this case 
taken as 4 /m ) and C-, T is required. For the original calculations, values of 
c„ _ rr were taken from the Monaghan T-prime method (appendix A and fig. 12) but any 

* jOj 1 

other reasonable values for classical flat-plate laminar skin friction could be used such 
as those, for example, from the Crocco method. In the present case to obtain Cp with- 
out viscous interaction for the triangular planform case, Cp Q rp must be multiplied by 
4/3. (See statements after eq. (43).) 

LIMITATIONS OF THEORY 

There are two regions on a plate where deviations from the local similarity theory 
might be expected. The first is very close to the leading edge and the second is the 
trailing-edge region. In the near leading-edge region the effect is noted as a departure 
from the viscous interaction theory with a tendency toward a pressure plateau which has 
been found in several experiments. There have been several interpretations of this effect 
but, in general, the parameter determining its onset is much the same from various theo- 
ries and appears to agree with experiment. Talbot (ref. 6) interprets this region of 
departure as an effect of slip flow and his result for the upstream limit of applicability of 
ordinary viscous interaction can be approximated as 
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(50a) 


*lim 


0.043 


y + 1 
y - 1 



or in terms of X 

X lim * 0.0185(y + 1) + 0.39J M J (50b) 

This result agrees with the limitation set by Oguchi (as given by Mann and Bradley, ref. 7) 
from different considerations. A graph of equation (50b) is presented as figure 17 applica- 
ble for y = 7/5; however, since the limit is not precise, the curves could be used for 
values of y from 1.2 to 1.67. Since aerodynamic forces are an integration of local 
values, the limits shown in figure 17 are considered to be optimistic for these forces. 

The second limitation is due to the upstream influence of the pressure changes at the 
wing trailing edge. This effect tends to reduce the aerodynamic forces. No definite 
limits have been set for this influence which is considered by Rogers and Metcalf in some 
detail in reference 8 for wings in supersonic flow. The complete pressure distributions 
from flat and curved plates in references 9 and 10 tested at Mach number 6.9 indicate a 
small influence of the trailing edge on the lower surface but a significant effect on the 
upper surface. The Mach numbers and Reynolds numbers for the tests of references 8, 

9, and 10 are widely different. 


CONCLUDING REMARKS 

Solutions have been obtained of the effect of viscous interaction not only on the pres- 
sure and skin friction but also on the aerodynamic forces of a sharp-leading-edge flat sur- 
face. The analysis assumes the applicability of laminar-hypersonic-local-similarity- 
boundary-layer theory in a perfect gas with Prandtl number unity and a constant ratio of 
specific heats. The results are presented in the form of correlation graphs which cover a 
large range of plate incidence and viscous interaction parameter and allow the effects of 
viscous interaction to be readily determined. Computations are shown for ratios of spe- 
cific heats of 7/5 and 5/3 and, for the aerodynamic forces, both the two-dimensional 
wing and the triangular planform wing are considered. The results indicate that the 
effect of viscous interaction can reduce the lift -drag ratio significantly. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 3, 1966. 


21 



APPENDIX A 


METHOD OF SOLUTION 


To solve the system of equations (1), (2) (or (3)), and (4) for obtaining the pressure 
distribution, the Runge-Kutta method to fourth order was used starting at the leading edge 
(X = °°) with the value predicted by strong interaction theory (eq. (11)) and progressing in 
the x-direction. The increments in X set into the high-speed digital computer were as 
follows: 


X 

CO 

VII 

< £ 

VII 

o 

AX 

(K 0 > 1*1) 

>100 

10 

5 

100 to 10 

5 

2.5 

10 to 1 

.025 

.00125 

1 to 0.1 

.0025 

.00125 

0.1 to 0.01 

.00125 

.000625 

0.01 to 0.001 

.000325 

.000125 


Because of the extreme sensitivity of pressure to small errors as the smallest values of 
X were approached, double precision (16 significant figures) was utilized for all phases 
of the calculation. Even this accuracy was not sufficient to obtain a solution at the smaller 
values of X. In these cases a second-order equation was fitted between the known values 
of dP/dX at X = 0 (eq. (5) or (6)) and the last value of dP/dX ^at X^j considered to 
be given accurately by the solution to the complete equations, and the solution was then 
carried out as before, except for the fixed values of dP/dX. The fit to dP/dX in these 
cases was 


dP 

dX 


= 3 


Ls - 2d N 

\ X F 3 ) 



+ 2 


/3d - 


x F g' 


IX + 


dP\ 

dx / X =0 


(Al) 


where 


d = P F - P 0 


+ x F 



g = 




The values of Xj, (the switching point) used with equation (48) were as follows: 
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K q = 15 

x F = 100 

8 S Kq S 12 

ii 

o 

CD 

VII 

VI! 

Xp - 1 

-0.75 SKg i3 

*1 

1! 

O 

CO 

Ml 

All 

rH 

1 

x F = 0.01 

CO 

V 

No switch 


As stated previously where the first portion of an integral was solved in closed form in 
order to provide an accurate start to a solution, the interval to which the closed form 
applied extended from X = «> to X = Xj = 10^ (eqs. (10), (13), (26), (39), and (43)). 

Where coefficients and constants were needed for solving for the aerodynamic coef- 
ficients, they were taken as follows: G which relates X and x is given by Monaghan 
(ref. 11) as 


1 /T 4.65Nt3 

G = 1.7208^-^— + — 

2 l TV 


1/3 


3.65N 


Pr 


1/2 


2.59 


■N- 


Pr 


1/2 


(A2) 


For a Prandtl number of unity, the result from equation (A2) agrees with that obtained 
from laminar-hypersonic- local-similarity-boundary-layer theory (ref. 1) where 


G 


1.7208 



— + 0.38591 


(A3) 


Equation (A3) is that used for the computations in this paper. 


The coefficient C in the linear formula for viscosity is according to the T-prime 
method 


C 


= J£L1~ 

M=o T' 


(A4) 


where the prime designates evaluation of the parameter at the reference temperature 
(T-prime). For the temperature ratio T' /Too Monaghan's constants were used so that 
(for N pr = 1) 

T 

= 0.468 + 0.532— + 0.195^=-^ M^ 2 (A5) 

Too Too 2 00 


Calculations were done for air and helium. For air, Keyes' three -constant modified 
Sutherland formula was assumed to apply for viscosity (ref. 12) and for helium a power 
law was assumed to represent the viscosity so that C was given by 
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APPENDIX B 

APPROXIMATE SOLUTIONS FOR VISCOUS INTERACTION 
PRESSURE DISTRIBUTION 


To first order for values of P greater than but not too close to 1, equation (2) may 
be written as 


P = 1 + r(r + D . K 2 

2 


(Bl) 


Combining equation (Bl) with equation (1) and with K = K 0 + Kg and K^ = 


= 1 


^prrJlIr^l 




/v. 

Ip 


(i + _*_dP\l 

\ + 2P dX/_ 


(B2) 


If furthermore dP/dX is assumed to be essentially constant so that dP/dX may be 
taken as (P - P Q )/X without large error, equation (B2) becomes 




/ 2P 

Vr(r + 1) 

3 -Lo 


where equation (B3) is applicable to both positive and negative Kq. 
If equation (B3) is restricted to Kq = 0, 


X = 2 


/ 2P _/ 


f Po - 1 \ 

Vr(y + D\ 

3-^2 / 


(B3) 


(B4) 


Equations (B3) and (B4) have the disadvantage that P is not an explicit function of 
X; however, they do show the proper behavior of P with X and equation (B4) exhibits 
the peak in (P - P Q )/X given by the "exact" solution (fig. 1). The accuracy, compared 
with the "exact" solution, increases as Kq increases. Although K^ has been assumed 
as unity in this solution, the value of (P - P Q )/x for a value of K^ other than unity may 
be obtained by simply multiplying the value of (P - P Q )/X for K^ = 1 by the proper 
value of Kj. 

Further approximations of equation (B4) allow the pressure parameters to be made 
a function of X. For one solution let P — P Q thus, 

= yfe r(y + l) (B5) 


2 ^ 


APPENDIX B 


Equation (B5) is a slightly less restrictive form of the "exact” result for Kq = <*> with 
A = 0 given after equation (5). It may be noted that the equation for Kq = 00 is equiva- 
lent to that given in reference 4. Equation (B5) is reasonably close to the exact result 
from equation (5) down to Kq ~ 0.5. 


Another solution may be obtained for P » P Q (large A), and equation (B4) 
becomes 



Equation (B7) is the same as equation (11), except for which can be put in as a multi- 

plication factor as discussed after equation (B4). 
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TABLE I.- VALUES OF SEVERAL USEFUL ASYMPTOTIC PRESSURE PARAMETERS 



\fPo 

y = 7/5 

oo 

19.498 

15.623 

13.045 

11.758 

10.472 

7.9136 

6.6435 

5.3852 

4.1483 

3.5440 

2.9553 

2.3907 

1.8635 

1.6198 

1.3931 

1.1859 

1.1089 

1.0718 

1.0354 

1 

.96544 

.93173 

.89888 

.83567 

.69159 

.56620 

.45795 

.28697 

.16731 

.088388 

.0404077 

.003578 

0 



2.5923 

2.5849 

2.5809 

2.5759 

2.5722 

2.5670 

2.5484 

2.5305 

2.4997 

2.4409 

2.3903 

2.3136 

2.1945 

2.0096 

1.8852 

1.7388 

1.5741 

1.5050 

1.4701 

1.4350 

1.4000 

1.3653 

1.3310 

1.2974 

1.2315 

1.0758 

.93256 

.80141 

.57395 

.39040 

.24749 

.14167 

.02504 

0 


Po 

OO 

502.25 

322.25 
224.47 
182.24 
144.46 

82.231 

57.779 

37.764 

22.180 

16.043 

11.000 

7.0355 

4.1142 

3.0225 

2.1562 

1.4919 

1.2762 

1.1781 

1.0862 

1 

.91940 

.84408 

.77378 

.64723 

.40188 

.23730 

.13169 

.031250 

.0041152 

.0001286 

0 


\/Po 
y = 5/3 

oo 

22.411 

17.951 

14.982 

13.500 

12.019 

9.0681 

7.6012 

6.1452 

4.7096 

4.0054 

3.3166 

2.6525 

2.0283 

1.7385 

1.4684 

1.2214 

1.1297 

1.0854 

1.0422 

1 

.95885 

.91874 

.87965 

.80451 

.63394 

.48713 

.36289 

.17678 

.064150 

.01134 

0 



2.9814 

2.9748 

2.9711 

2.9666 

2.9632 

2.9584 

2.9414 

2.9249 

2.8960 

2.8399 

2.7905 

2.7136 

2.5896 

2.3870 

2.2448 

2.0727 

1.8756 

1.7924 

1.7504 

1.7085 

1.6667 

1.6252 

1.5840 

1.5432 

1.4627 

1.2679 

1.0825 

.90722 

.58926 

.32075 

.11340 

0 
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(a) r = 7/5; windward surface. 


Figure 3.- Average pressure parameter as a function of viscous interaction parameter for a two-dimensional 
flat plate at various inclinations to undisturbed flow. 
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Figure 7.- Effect of variations in the viscous interaction parameter on initial normal-force-curve slope. Zero thickness plate. 
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(a) r = 7/5; T w /T t = 1. 

Figure 10.- Effect of viscous interaction on average skin-friction contribution of the windward and leeward surfaces 

of a two-dimensional fiat plate. 














(e) r = 7/5; T w /T t = 0 . 2 . 
Figure 10 .- Continued. 





































(g) r = 5/3; T w /T t = 1 . 
Figure 11.- Continued. 
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(hi r = 5/3; T w At = 0.8. 
Figure 11.- Continued. 
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(I) r = 5/3; T w /T t = 0 . 
Figure 11.- Concluded. 
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(a) Air. 

Figure 12.- Total average skin friction without viscous interaction for a two-dimensional flat plate at zero angle of attack utilizing 
Monaghan's method. Sharp leading edge and Prandtl number 1. 
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(a) y = 7/5. 

Figure 16.- Effect of viscous interaction on ratio of normal-force-curve slope to average skin friction at 
zero angle of attack for zero-thickness plates. 
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Figure 17.- Upstream limit for applicability of local similarity solution. Curves shown are for 
Y = 7/5 but for values of y between 1.2 and 5/3 they change only about 10 percent. 
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